Dark radiation and dark matter coupled to holographic Ricci dark energy 
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We investigate a universe filled with interacting dark matter, holographic dark energy, and dark 
radiation for the spatially flat Friedmann-Robertson- Walker (FRW) spacetime. We use a linear 
interaction to reconstruct all the component energy densities in terms of the scale factor by directly 
solving the balance's equations along with the source equation. We apply the x 2 method to the 
observational Hubble data for constraining the cosmic parameters, contrast with the Union 2 sample 
of supernovae, and analyze the amount of dark energy in the radiation era. It turns out that our 
model exhibits an excess of dark energy in the recombination era whereas the stringent bound 
fi x (z ~ 10 10 ) < 0.21 at big-bang nucleosynthesis is fulfilled. We find that the interaction provides 
a physical mechanism for alleviating the triple cosmic coincidence and this leads to $! m o/fi x o — 
fi r o/fi x o - 0(1). 
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I. INTRODUCTION 

The two main components of the universe are dark 
matter and dark energy. The dark matter accounts for 
1/3 of the stuff in the universe and is also inextricably 
connected with the formation of galaxies and galaxy clus- 
ters. In fact this new form of matter not only holds galax- 
ies together, but also is responsible for the large-structure 
formation in the universe [l| . The astrophysical evidence 
for dark matter come form colliding galaxies, gravita- 
tional lcnsing of mass distribution or power spectrum of 
clustered matter Q. 

The other 2/3 exists as in an even more mysterious 
form dubbed as dark energy and is causing the expan- 
sion of the Universe to speed up, rather than slow down. 
The first serious observational hints of dark energy in the 
universe date back to the 1990s when astronomers obser- 
vations of supernova were used to trace the expansion 
history of the universe [|| . Subsequently, the WMAP re- 
sults suggested that the aforesaid amount of dark energy 
could explain both the flatness of the universe along with 
the observed accelerated expansion 0]. Nowadays, there 
is a growing number of observational methods for prob- 
ing the dynamical behavior of dark energy at different 
scales; galaxy redshift surveys allow to obtain the Hub- 
ble expansion history by measurement of baryon acoustic 
oscillation in the galaxy distribution [l[, geometric weak 
lcnsing method applied to Hubble space telescope images 
helps to find tighter constraints on the dark energy equa- 
tion of state also [f|. 

Despite the devoted effort for understanding the nature 
of dark matter and dark energy, there was not found a 
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microscopic theory for the dark side of the universe, ca- 
pable of unraveling their particles content. On the ob- 
servational side, the data of several tests have confirmed 
some of their plausible properties, such as that it could 
be the repulsive effect of dark energy or the clustering 
action of dark matter. 

Another interesting trait to explore is related with 
the existence of non-gravitational coupling between dark 
matter and dark energy, such exchange of energy could 
alter the cosmic history, leaving testable imprints in the 
universe Jj|. It is believed that a coupling between dark 
energy and dark matter changes the background evolu- 
tion of the dark sector allowing to constrain any type 
of interaction and giving rise to a richer cosmological dy- 
namics compared with non- interacting models @ . A step 
forward for constraining dark matter and dark energy 
is to use the physics behind recombination or big-bang 
nucleosynthesis epochs by adding a decoupled radiation 
term to the dark sector for taking into account the strin- 
gent bounds related to the behavior of dark energy at 
early times Q, @|- The behavior of dark energy in the 
recombination era was explored within the framework of 
three interacting components also [9|. 

Our goal is to consider a model where dark matter 
and dark radiation are coupled to holographic dark en- 
ergy and explore the cosmic triple coincidence problem 
related to the amount of these components at present 
fioj . We perform a cosmological constraint using the up- 
dated Hubble data [ll|, numerically obtain the distance 
modulus fi(z) for contrasting with the Union 2 compi- 
lation of supernovae la and analyze the order of 
magnitude corresponding to the cosmological parameter 
known as transition redshift. In order to check the fea- 
sibility of the model, we also examine the severe bounds 
for dark energy in the recombination era (lH ] or nucle- 
osynthesis epoch (l4T |. 
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II. THE MODEL 

We consider a spatially flat homogeneous and isotropic 
universe described by FRW spacctimc with line element 
given by ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) being a(t) 
the scale factor. The universe is filled with three inter- 
acting fluids namely, dark radiation, dark matter and 
modified holographic Ricci dark energy so that the evo- 
lution of the FRW universe is governed by the Friedmann 
and conservation equations, respectively, 



'SH 2 = p T + Pn 



(1) 



p + 3H(p T + p r + p m + p m + p x + p x ) = 0, (2) 

where a is the scale factor and H = a/a stands for the 
Hubble expansion rate. Here, we will use the holographic 
principle within the cosmological context by associating 
the infrared cutoff L with the dark energy density, thus 
we take L~ 2 in the form of a linear combination of H and 
H 2 : 



H+'^H 2 
P \ 2 



(3) 



Here, a and j3 are two free constants. In particular, we 
obtain p x oc R for a = 4/3, where R = Q(H + 2H 2 ) is 
the Ricci scalar curvature for a spatially flat FRW space- 
time. 

The use of the variable t] = ln(a/ao) 3 , where ao is set 
as the value of the scale factor at present and ' = d/drj, 
allows us to rewrite Eqs. ([2]) and (|3|) as 



P = ~7rPr 



7m Pn 



7xPx, 



P = -upc - ftp*, 



Pc = Pr + Pn 



(4) 



(5) 



(0) 



where ji = I + Pi/pi denotes the barotropic index, 
not necessarily constant, of each component with i = 
{r, m,x}, 7 r ~ 4/3, 7, n ~ I, < 7 X < 2/3 so that 
< 7x < 7m < 7r- Taking into account Eq. (|5|) along 
with Eq. ©, we could extract as a physical hint that 
the modified holographic dark energy ([3]) forces to the 
dark matter and dark radiation to have the same bare 
equation of state. 

The holographic dark energy © or ([5]), looks like a 
"conservation equation" for the three dark components 
with constant coefficients. Therefore, the selected holo- 
graphic dark energy ([3]) or ((5j) has transformed the orig- 
inal model of three interacting components into another 
simpler scheme of two interacting components having two 
constant equations of state . As we have already men- 
tioned above, such degeneration occurs because the dark 
radiation and dark matter have the same equation of 
state, p r = {a — l)p r and p m = (a — l)p m ■ After com- 
paring the whole conservation equation Q with modified 



conservation equation (|5|), we obtain the compatibility 
relation 



li-Pr + ImPm + 7xPx = a(p r + p m ) + j3p x , 



(7) 



that relates the equation of state of the dark components 
with the bare ones. In what follows, we will use Eq. 
([5]) with constant coefficients a and /3 instead of Eq. ^ 
with the non-constant coefficient 7 X = l+p x /p x . In some 
sense, Q and ((SJ) give rise to different representations of 
the mixture of two interacting fluids and clearly these de- 
scriptions are related between them by the compatibility 
relation ([?])• Therefore, the holographic dark energy Q 
conveniently links a model of three interacting fluids hav- 
ing variable equations of state with a model of two 
interacting fluids with "bare constant equations of 
state" . 

Solving the system of equations p = p c + p x and ([5]) we 
get the energy densities of both component as a function 
of p and its derivative p' 



0P + P* 



ap + p 



(8) 



A A ' 

where A = a — f3 is the determinant of the linear equation 
system and we assume that j3 < a. Now, we introduce 
energy transfer between those components by splitting 
Eq. ([5]) into two balance equations 



Pc + "Pc 



PL + PP* = Q, 



(9) 
(10) 



where we have considered a coupling with a factorized 
H dependence 37JQ, being Q the interaction term that 
generates the energy transfer between the two fluids. Af- 
ter differentiating the first Eq. ([8]) and combining with 
Eq. ([5]), we obtain a second order differential equation 
for the total energy density (source equation) 



p" + (a + f3)p' + af3p = QA. 



(11) 



We will take into account interactions Q for which the 
solutions of the evolution equation for the scale factor 
H = \fp~Ji includes power law ones, because they play 
an essential role for determining the asymptotic behav- 
ior of the effective barotropic index 7 = (ap c + /3p x )/p = 
—2H/3H 2 . It describes a universe approaching to a sta- 
tionary stage 7 S associated with the constant solution 
7 = 7s, where < (3 < 7 S < a and a = t 2 ^ ls . In addi- 
tion, given a set of initial conditions, if 7 tends asymp- 
totically to the constant solution 7 S , then 7^ becomes 
an attractor solution [f|. An interaction satisfying this 
requirement belongs to the class 



(7 S - op (7 S - 0) 
Q = 7 P, 



(12) 



with Q < . Solving the source equation (fTTj) for this 
interaction, we obtain the total energy density 



ca 



" 3 T' +ba-^ a+ ^\ 



(13) 
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Hence, for any initial conditions c, 6, and large scale fac- 
tor, the energy density behaves as c/a 37s and the power- 
law expansion a — > t 2 / 3ls becomes asymptotically stable. 

The dark densities of coupled components, p x and p c , 
are given by 

(7, - g) ca~ 3 ^ + (a- 7s )ba~ 3( - a +P-^ 

Pc = 7 , (14) 



Px = ^ ■ (15) 

Thus, the ratio r = p c /p x tends to r s = (7 S — j3)/(a — j s ), 
being r s an attractor. 

With the aid of the total energy density (fl"3)) , we can 
calculate the explicit form of the interaction term (fT2|) as 
a function of the scale factor, 



(is - a) (7 S - fi) 



ca-^'+ba- 3 ^-^ 



In turn Eq. ^ can be rewritten as 

Pr + Pm + a (Pr + Pm) ' 



(16) 



(17) 



In order to break the degeneracy of this set of compo- 
nents, formed by dark radiation and dark matter, we 
introduce partial interactions into the corresponding bal- 
ance equation of both components as follows: 



Pr + ap r 



apm = Qr, 



(18) 



(19) 



where Q r and Q m stand for the exchange of energy be- 
tween p r and p x , and besides these satisfy the condition 



+ < 



Qr =0, 



(20) 



to recover the conservation equation ([5]) after having 
summed all Eqs. (fT7} -(JT9" |) . Also we assume that Q r and 
Q m are a linear combination of the two terms contained 
in the interaction term Q (|16[) , so they read 



Q r = ci a 37s + &i a 



_3( Q+(3 _ 7s ) 



(21) 



Qm = -(ci + Q c) a- 3 ^ - (h + Q b) a~ 3 ^ + ^\ (22) 

where Qo = ( 7s — a )(js — /3)/A, while c\ and b\ are 
free parameters of the model. Inserting these interac- 
tions into the evolution equation of the dark radiation 
and dark matter (|18p - (fTl?)) and solving these coupled sys- 
tem of equations we obtain 



Pr = 



aa- 3 ^ b ia - 3 ^ a+f} -^^ 



a - Is 



— Co a 



Is -13 



(23) 



( Cl +Q c)a~ 3 ^ _ 3a ( Ol + Q &) a - 3 («+^) 
-+cq a 



a - 7 S 



(24) 



where the first and third terms in both dark energies den- 
sities are the particular solutions of the evolution equa- 
tions (fT8f and (fl~9|) while the second term, in both Eqs. 
(f2"3"|) and (fM]) . is the homogeneous solution of the linear 
system of equations ([18]) and (fl~9)). In what follows, we fix 
Ci = — Qoc/2 > and &i = — Qob/2 > without loss of 
generality. Using Eq. (|T3j). (HHJ), (ESJ), and ([M]), we find 
the coefficients c, 6, and cq in terms of the the density 
parameters and 7 S . Now, we only show their expressions 
for a = 4/3 and a + (3 — 7s = 1: 



c = 3iJo 
6 = ZH'l 
c = 3ff 2 



l + n x0 (3 7s -5) 



3(7s - 1) 
(50 x o-4) + 3 7s (l-n x0 ) 



3(7, - 1) 
1 — ^xo — 2f2 m o 



(25) 
(26) 
(27) 



Here, fiio stands for the density parameter of each com- 
ponent. 



III. OBSERVATIONAL CONSTRAINTS ON THE 
THREE INTERACTING MODEL 

We will provide a qualitative estimation of the cosmo- 
logical parameters by constraining them with the Hub- 
ble data [l5|- [l(| and the strict bounds for the behavior 
of dark energy at early times [lj]. In the former case, 
the statistical analysis is based on the x 2 -function of the 
Hubble data which is constructed as (e.g. |T3|) 



x 2 w = E 



[H(9,z k )-H ohs {zk)Y 



k=l 



a(zkf 



(28) 



where 9 stands for cosmological parameters, i/ob s (zfc) is 
the observational H(z) data at the redshift cr(zk) is 
the corresponding lcr uncertainty, and the summation is 
over the 19 observational H(z) data fill]. Using the abso- 
lute ages of passively evolving galaxies observed at differ- 
ent redshifts, one obtains the differential ages dz/dt and 
the function H (z) can be measured through the relation 
H(z) = -{l + z)- 1 dz/dt [H, [3. The data H ohs ( Zi ) and 
H i, s (zk) are uncorrelated because they were obtained 
from the observations of galaxies at different redshifts. 

From Eq. (fT3"f one finds that the Hubble expansion of 
the model becomes 



H(6\z) = H (Cx 3 ^ + Bx 3( > a+ P- 



7a) 



(29) 



where c — 3HqC, b = 3HqB as obtained form (|25|) . 
(f2"6"|) .respectively. Here, we consider 9 = {i7 ,7s,^xo} 
as the independent parameters to be constrained for the 
model encoded in the Hubble function (|29|) with the sta- 
tistical estimator (|28l) . while a is taken equal to 4/3 to 
have an early era dominated by radiation. Besides, we 
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2D Confidence level 



Priors 


Best fits 


Xd.o.f 


7s = 10" 3 
tt x0 = 0.67 
Hq = 70.4 


(ifo,n x0 ) = (72.05±! 9 9 », 0.637™) 
(ff o , 7s ) = (73.8li^t0.00lt°;*°t) 
(n x0l7s ) = (0.623™, 0.001™) 


0.81 
0.94 
0.86 



TABLE I: Observational bounds for the 2D C.L. obtained in Fig. 
( TT1 by varying two cosmological parameters. The Xd o f m a ^ ^ ne cases 
studied is less the unity; in fact it goes from 0.81 to 0.94. 

will impose a + j3 — 7 S = 1 so that the universe ex- 
hibits an intermediate stage dominated by pressureless 
dark matter. For a given pair of (0\, #2), we are going to 
perform the statistic analysis by minimizing the x 2 func- 
tion to obtain the best fit values of the random variables 
#crit = {^criti) $crit2} which correspond to a minimum of 
Eq. (P5|) . Then, the best-fit parameters 6* cr it are those val- 
ues where Xmm(^crit) leads to the local minimum of the 
X 2 {6) distribution. If xl.o.i = X„ lin (^rit)/(A r - n) < 1 
the fit is good and the data are consistent with the con- 
sidered model H(z;8). Here, N is the number of data 
and n is the number of parameters [I?} ■ The variable % 2 
is a random variable that depends on N and its proba- 
bility distribution is a x 2 distribution for N — n degrees 
of freedom. Here N = 19 and n = 3, so in principle, 
we will perform three minimizations of the x 2 statistical 
estimator, interpreting the goodness of fit by checking 
the condition Xd f < 1> as a WS V to keep the things 
clear and focus on extracting relevant physical informa- 
tion from this statistical estimation. 
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FIG. 1: Two-dimensional C.L. associated with ler,2cr for different 6 
planes. 



The random datasets that satisfy the inequality Ax 2 = 
X 2 ($) — Xmin(^crit) < 2.30 also represent 68.3% confidence 




FIG. 2: (Left panel) :Thrcc-dimcnsional C.L. for the Hq — 7 S — f7. x o 
plane. (Right panel): Two-dimensional C.L. obtained after have per- 
formed the marginalization over Hq. 
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FIG. 3: Hubble diagram for the Union2 compilation: the points 
are the observational data of supernovae, the theoretical curves fi(z) 
(dashed lines) represent the best cosmological model for different 
cases: (global estimation) / - (H , fi x0 , j s ) = (72.05,0.63,0.001), 
(marginalizing over H a ) II - (O x0 ,7 s ) = (0.6191° °,H , O.OOll °' 2 2 H) , 
III - (ffo,n x0 ) = (72.05+2 l° 3 , 0.6371°,'?!!). Iv ~ (#0,7=) = 
(73.8lil 2 t,0.00lir 9 ° 8 t), V - (O x o,7 3 ) = (0.623™ • 0-001™)- 



contours in the 2D plane at lcr level. It can be shown 
that 95.4% confidence contours with a 2cr error bar in 
the samples satisfy A% 2 < 6.17. The two-dimensional 
C.L. obtained with the standard x 2 function for two in- 
dependent parameters is shown in Fig. ([1]), whereas the 
estimation of these cosmic parameters is briefly summa- 
rized in Table dJ). 

We obtain j s < 1CP 3 , so these values clearly ful- 
fill the constraint 7 S < 2/3 which ensures the exis- 
tence of an accelerated phase of the universe at late 
times [Table ©]. We find the best fit at (H ,n x0 ) = 
(72.05t 2 ; 9 9 °km s" 1 Mp C -\ 0.6371 ,; ,!) with xL.f = 0-81 
by using the prior j g = 10~ 3 . These findings show, in 
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FIG. 4: (Upper panel) :Thrcc-dimonsional C.L. for the Ho — —7s — 
— z c plane. (Lower panel): Two-dimensional C.L. obtained after have 
performed the marginalization over Hq. 
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FIG. 5: Plots of Q x (z), Q m (z), O r (z), q(z), uj c f{(z), uj e ff x (z) using 
the best-fit values obtained with the Hubble data for different 9 planes. 
Plot of LogO x in terms of Log z. 



broad terms, that the estimated values of Ho and f2 xu 
are in agreement with the standard ones reported by the 
WMAP-7 project [ljj. The value of Sl x0 is slightly lower 
than the standard one of 0.7 being such discrepancy less 
or equal to 0.1%. We find that using the priors Ho = 
70.4km s" 1 Mpc" 1 the best-fit values for the present- 
day density parameters are (S7 x o , £l r o) = (0.62, 7 x 10~ 2 ) 
along with a larger goodness condition [x\ f = 0-86) 
[Tabic (jj)]. Regarding the amount of dark matter at 
present, we have fixed f2 m o = 0.3 because this value is 
consistent with one reported by the WMAP-7 project 
[l9l ]. In performing the statistical analysis, we find that 
H 6 [70.4, 73.81]km s" 1 Mpc -1 so the estimated values 
arc met within la C.L. reported by Ricss et al [Q, to 
wit, H Q = (72.2 ± 3.6)km s _1 Mpc _1 . In order to en- 
sure that our previous local minimization analysis is cor- 



rect, we have performed a global statistical analysis by 
estimating all the parameters at once. In doing that, 
we obtain (H , ft x0 , 7s) = (72.05, 0.63, 10^ 3 ) along with 
x l o f = 0.86 < 1 [see Fig. ©]. Fig. © shows two- 
dimensional C.L in the H x o 7 S plane obtained when 

the joint probability P(Hq,Q x q,^/ s ) is marginalized over 
Ho (see Fig. ©); then the marginalized best-fit val- 
ues become (ft x0 ,7 s ) = (0-619±g;g||, 0.001±g;H|) with 
Xd o f = < 1- ^ must be stressed that we report 
for the most relevant minimization procedures the corre- 
sponding marginal la error bars [201 ] as can be seen in 
Table ©. 

As is well known, distance indicators can be used for 
confronting distance measurements to the corresponding 
model predictions. Among the most useful ones are those 
objects of known intrinsic luminosity such as standard 
candles, so that the corresponding comoving distance can 
be determined. That way, it is possible to reconstruct the 
Hubble expansion rate by searching this sort of object 
at different redshifts. The most important class of such 
indicators is type la supernovae. Then, we would like 
to compare the Hubble data with the Union 2 compila- 
tion of 557 SNe la [l2j by contrasting theoretical distance 
modulus with the observational dataset. In order to do 
that, we note that the apparent magnitude of a supernova 
placed at a given redshift z is related to the expansion 
history of the Universe through the distance modulus 



M = 5 log 



d L {z) 



Mo, 



(30) 



where m and M are the apparent and absolute magni- 
tudes, respectively, (io = 42.38, h = iJo/100km _1 s~ 1 , 
and cLl{z) = H {1 + z)r{z), r(z) being the comoving dis- 
tance, given for a FRW metric by 



r(z) 



dz' 

W) 



(31) 



Using the Union 2 dataset, we will obtain five Hubble 
diagrams and compare each of them with the theoreti- 
cal distance modulus curves that represent the best-fit 
cosmological models found with the update Hubble data 
(see Fig. ©); it turned out that at low redshift (z < 1.4) 
there is an excellent agreement between the theoretical 
model and the observational data. 

For the sake of completeness, we also report bounds for 
other cosmological relevant parameters [see Table ([H])]. 
such as the fraction of dark radiation f2 r (z = 0), the 
effective equation of state at z = (w e ffo = 7 e gfo — 1), de- 
celerating parameter at the present time qo, and the tran- 
sition redshift (z c ) among many others, all these quanti- 
ties are derived using the three best fit values reported in 
Table (UJ) . We find that the z c is of the order unity vary- 
ing over the interval [0.68, 0.931 , such values are close to 
z c = 0.69^0 13 reported in [2l|, [22j quite recently. More- 
over, taking into account a x 2 -statistical analysis made in 
the (u!q, z c )-plane based on the supernova sample (Union 
2) it has been shown that at 2a C.L. the transition red- 
shift varies from 0.60 to 1.18 23[. In order to estimate 
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Bounds for cosmological parameters 


u 


:., 


9(2 = 0) 




-HI) 


u>„4z = o) 


ft,(z ~ 1100) 


f!,(z ~ 10 10 ) 


<l 


n, 


f2,„o/ft,o 


Si,o/S!*n 


1 


0.75 


-0.59 


-0.64 




-1.32 


0.26 


0.26 


0.3 


0.062 


0.47 


0.09 


II 


0.93 


-0.67 


-0.74 




-1.23 


0.23 


0.23 


0.3 


0.03 


0.44 


0.04 


III 


lies 


-0.55 


-0.70 




-1.28 


0.20 


0.20 


0.3 


0.07 


0.48 


0.1 



TARLE II: Derived bounds for cosmic parameters using the best fits 
value of 2D C.L. obtained in Table. Q by varying two cosmological 
parameters in three different cases. 



z c as independent parameter using a % 2 method, we first 
needed to obtain its generic formula by imposing the con- 
dition (z c + l)H'(z c ) = H(z c ): 



30,0(1 - 7,) + (4 - 3 7s ) 
(2-37 fl )[3J2 x o(l-7«) + l] 



1/3(1-7.0 



1. (32) 



Hence, placing Eq. ([32]) into Eq. ([29]), the Hubble 
parameter turns out to be a function of Hq, z c , and 
7 S . The global statistical analysis using z c as indepen- 
dent parameter instead of fi^o leads to (Hq, z c ,~/ s ) = 
(72.05, 0.63, 10" 3 ) along with xL.f = °- 86 [ see Fig- 
ill])] whereas the marginalized best-fit values become 
(z c , ls ) = (0.623t° : °52 5 0.00lt°; 3 7 ^) with a xL.f = 0.86 
[see Fig. Q]; notice that the marginalized best fit value 
of z c is considerably lower than the values reported in 
Table pi]). Besides, the behavior of decelerating param- 
eter with redshift is shown in Fig. ([5]), in particular, its 
present-day value varies as —0.67 < qo < —0.55 for the 
three cases mentioned in Table II, and all these values are 
in perfect agreement with the one reported by WMAP-7 
project fl~9| . 

The effective EOS of the mix is given by 



W ff 



/3f2 x + Q.'il m + C*^r 



i. 



(33) 



The effective equation of state (EOS) for dark energy is 
obtained from Eq. (fTTJj) . it reads 



Weffx = 7x 1- 

Px 



(34) 



In Fig. ([5]) we plot the effective equation of state as a 
function of redshift for the best-fit value shown in Table 
I, in general, we find that w er j > — 1 provided that (1 — 
37 s )f2 x + 4(O r + fl m ) < 3, as a matter of fact its present- 
day values cover the range [—0.74, —0.64]. On the other 
hand, the effective EOS associated to the dark energy 
evaluated at z = 0, lo s x (z = 0) varies over the range 
[-1.32,-1.23]. 



In regard to the behavior of density parameters fi x , 
f2 m , and f2 r , we see that very close to z = the dark 
energy is the main agent that speeds up the universe, 
far away from z = 1 the universe is dominated by the 
dark matter and at very early times the radiation compo- 
nent governs the entire dynamic of the universe around 
z ~ 10 3 [cf. Fig. © ]. In this point, we would like 



to present an appealing discussion concerning the triple 
cosmic coincidence problem (TCC) [T3| related to the 
amount of dark energy, dark matter, and radiation at 
present moment. We have proposed a physical mecha- 
nism based on a phenomenological interaction among the 
three cosmic components, in fact, since we are working 
within the framework of three interacting cosmic compo- 
nents the aforesaid scenario seems to be a fertile arena for 
studying the TCC. As pointed out by Arkani-Hamed et 
al in their seminal work: " there is an era in the history of 
the universe where all three forms of energy, in matter, 
radiation and dark energy, become comparable within 
a few orders of magnitude" fl0| . Here, we have found 
that interaction made it possible to have f2 m o/^xO = 
0.3/0.62 = 0.48 and f2 r0 /f2 x0 = 0.07/0.62 = 0.11, so 
f2 m o/r2 x o — ^ro/^xo — C(l)i showing that the interac- 
tion implemented can be used for alleviating the TCC 
[see Table ©]. 

Now, we explore another kind of constraint which 
comes form the physics at early times because this can 
be considered as a complementary tool for testing our 
model. As is well known the fraction of dark energy in the 
recombination epoch should fulfill the bound f2 dc(z — 
1100) < 0.1. Taking into account the best-fit values re- 
ported in Table (PJ, we find that at early times the dark 
energy does not change much with the redshift z over the 
interval [10 3 , 10 10 ], in fact, the Log ft x in terms of Log z 
goes from 0.64 to 0.20 [sec Fig. ©)]. Table (HU shows 
that around z ~ 1100 (recombination) f2 x can vary from 
0.20 to 0.26. This excess of dark energy requires further 
research because some signal could arise from this early 
dark energy(EDE) models uncovering the nature of DE 
as well as their properties to high redshift, giving an in- 
valuable guide to thephysics behind the recent speed up 
of the universe [l3| , [|| ■ Regarding the values reached by 
fi x around the big bang nucleosynthesis (BBN) z = 10 10 , 
we find that there is a variation from 0.20 to 0.26, so the 
conventional BBN processes that occurred at tempera- 
ture of lMev is not spoiled because the severe bound 
reported for early dark energy f2 x (z ~ 10 10 ) < 0.21 is 
marginally fulfilled at BBN [l4|. 

As is well known, dark energy dominates the whole dy- 
namics of the universe at present and there is an obvious 
decoupling with radiation practically. However, from a 
theoretical point of view, it is reasonable to expect that 
dark components can interact with other fluids of the 
universe substantially in the very beginning of its evo- 
lution due to processes occurring in the early universe. 
For instance, dark energy interacting with neutrinos was 
investigated in [24j . The framework of many interact- 
ing components could provide a more natural arena for 
studying the stringent bounds of dark energy at recombi- 
nation epoch. There could be a signal in favor of having 
dark matter exchanging energy with dark energy while 
radiation is treated as a decoupled component 0|, or 
the case where dark matter, dark energy, and radiation 
exchange energy. More precisely, when the universe is 
filled with interacting dark sector plus a decoupled ra- 
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diation term, it was found that Q x (z a 1100) = 0.01 
or ft x (z ~ 1100) = 10~ 8 [| but if radiation is cou- 
pled to the dark sector, the amount of dark energy is 
drastically reduced, giving Q x (z ~ 1100) ~ ©(lO^ 11 ) [|. 
In our model, we have found that the amount of early 
dark energy varies in the range [0.26; 0.20], so the behav- 
ior of dark energy at recombination is considerably much 
smoother than in the aforesaid cases Q, 0, 0- We ex- 
pect to include a (decoupled) neutrino term in Friedmann 
equation to examine in more detail the dark radiation as 
a signature of dark energy. 

IV. CONCLUSION 

We have discussed a class of interacting dark matter, 
dark radiation, and holographic Ricci-like dark energy 
model for a spatially flat FRW background. We have 
coupled those components and obtained their energy den- 
sities in terms of the scale factor. 

We have examined the previous model by constrain- 
ing the cosmological parameters with the Hubble data 
and the well-known bounds for dark energy at recombi- 
nation era. In the case of two-dimensional (2D) C.L., 
we have made three statistical constraints with the Hub- 
ble function [see Fig. (JTJ) and Table ©]. We have 
found that 7 S < 10~ 3 , so these values fulfill the con- 
straint 7 S < 2/3 for getting an accelerated phase of 
the universe at late times. We find the best fit at 
(#o,n*o) = (72.05l2;^kms- 1 Mpc- 1 ,0.637lo;o38) with 
Xd o f = 0-81 by using the prior j s = 10~ 3 . It turned out 
that the estimated values of Hq and f2 x o are in agree- 
ment with the standard ones reported by the WMAP- 
7 project [l9j |. Besides, we have found that Hq G 
[70.4, 73.81]km s" 1 Mpc" 1 , so the estimated values are 
met within lcr C.L. reported by Riess et al [l8j], to 
wit, H = (72.2 ± 3.6)km s" 1 Mpc" 1 . After having 
marginalized the joint probability P(Hq, f2 x o, 7 S ) over Hq 
[see Fig. (JTJ], we saw that the marginalized best- fit 
values are (n x0 ,j s ) = (0.619±g;g|§, 0.001±g;i||) with a 
x \ o f = 0.87 < 1 [see Fig. ©]. Using the best fits men- 
tioned in Table (|T|), we have numerically obtained the 
distance modulus of the supernova as predicted by the 
theoretical model and compare with the Union 2 dataset, 
finding that at low redshift (z < 1.4) there is excellent 
agreement between the theoretical model and the obser- 
vational data [see Fig. ||3J]. 

Regarding the derived cosmological parameters, for in- 
stance, the transition redshift z c turned out to be of 
the order unity varying over the interval [0.68,0.93], 
such values are in agreement with z c = 0.69^1® re- 
ported in [2l[-[22| , and meets within the 2a C.L ob- 
tained with the supernovae (Union 2) data in [23|. We 
also have performed a global statistical analysis using z c 
as independent parameter instead of £I x q which lead to 
(flo,2c,7«) = (72.05,0.63, lO" 3 ) along with X & . Q .f = °- 86 
[see Fig. Q], whereas the marginalized best- fit values 



are (z C7 j s ) = (0.6231^1, 0.OOII^) together with a 
Xdof = 0-^6 [see Fig. (0J]. Besides, with the decelerat- 
ing parameters q(z = 0) G [—0.67, —0.55] for the three 
cases mentioned in Table pi)), all these values are per- 
fectly i n ag reement with the one reported by WMAP-7 
project [19| [see Fig. (|SJ)]. 

Concerning the effective equation of state, we have 
found that cj e g > —1 and its present-day values vary 
over the ranges [-0.74, -0.64] [see Table |n} and Fig. 
(O]. The equation of state associated with dark energy 
satisfies the inequality w ff x < — 1. 

Besides, we have found that the fraction of dark ra- 
diation at present moment fi r o varies in the interval 
[0.03,0.07] for the three cases mentioned in Table (|Hl) . 
The dark energy amount f2 x (z) governs the dynamic of 
the universe near z = 0, whereas far away from z = 1 
the universe is dominated by the fraction of dark mat- 
ter Q m (z) and at very early times the fraction of radia- 
tion Q r (z) controls the entire dynamic of the universe 
around z ~ 10 3 [cf. Fig. (|SJ|]. We also have exam- 
ined the triple cosmic coincidence problem fioj within 
the framework of three interacting cosmic components, 
finding that interaction used in this work provides a phc- 
nomenological mechanism for alleviating TCC, leading 
to Ono/Qxo - ^ro/^xo ^ 0(1) [see Table (TH])]. 

Finally, we have found that at early times the dark en- 
ergy does not change much with the redshift z over the 
interval [10 3 , 10 10 ], in fact, the Log ft x in terms of Log z 
goes from 0.64 to 0.20 [see Fig. ©)]. Table (|TTJ| shows 
that around z ~ 1100 (recombination) H x can vary from 
0.20 to 0.26. The latter results indicate an excess of dark 
energy so it requires further research [l3[, Q , in fact 
it could be related with the degeneracy presents in the 
equation of states of dark matter and dark radiation. In 
order to explore this issue in more detail, we expect to 
include an additional (decoupled) neutrino term in Fried- 
mann equation; thereby, we will seek to distinguish the 
radiation term coupled to dark matter, where both com- 
ponent share the same bare equation of state, from the 
decoupled neutrino term. However, it must be stressed 
that the values reached by i7 x around the big-bang nu- 
cleosynthesis (BBN) z = 10 10 vary from 0.20 to 0.26, 
so the conventional BBN process is not spoiled because 
our estimations, in most of the cases mentioned above, 
fulfill the severe bound reported for early dark energy: 

n x ( z ~ io 10 ) < o.2i [3. 
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